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1 Introduction 



Haldane-Shastry (HS) spin chain is a well known quantum integrable model, where 
equally spaced spins on a circle interact with each other through pairwise exchange 
interactions inversely proportional to the square of their chord distances. Study of 



fact that the exact ground state wavefunction of this model coincides with the C/ — > oo 
limit of Gutzwiller's variational wave function for the Hubbard model, and also with 
the one-dimensional version of the resonating valence bond state proposed by Anderson 
[1,2]. Remarkably, HS spin chain can be explicitly solved in much greater detail than 
integrable spin chains with short-range interactions, has a Yangian quantum group 
symmetry and interestingly shares many of the characteristics of an ideal gas, but with 
fractional statistics [3-5]. The Hamiltonian of SU{m) HS model with N number of 
lattice sites is given by 



where — ji^/N and Pj^ is the exchange operator interchanging the 'spins' (taking m 
possible values) on j-th and k-th lattice sites. 

By using the motif representations associated with y (giZ^) Yangian symmetry of HS 
spin chain (1.1), one can find out its complete spectrum including the degeneracy factor 
for each energy level [6-8]. However, in practice, the computation of such degeneracy 
factors becomes very cumbersome for m > 2 and large values of N. Therefore, it is 
difficult to express the partition function of HS spin chain in a simple form (for arbitrary 
values of N and m) with the help of motif representations. Due to this reason, it is 
worthwhile to explore other approaches for calculating the partition function of HS spin 
chain. In fact, a rather simple expression for the exact partition function of SU{m) HS 
spin chain has been obtained recently [9] by applying the so called freezing trick [10- 
12]. This freezing trick utihzes the connection between SU{m) HS spin Hamiltonian 
and SU{m) spin Calogero-Sutherland (CS) model which has dynamical as well as 
spin degrees of freedom. More precisely, one takes the strong coupling limit of spin 
CS Hamiltonian, so that the particles freeze at their classical equilibrium positions of 
the scalar part of the potential and spins get decoupled from the dynamical degree 
of freedom. As a result, one can derive the partition function of HS spin chain by 
'modding out' the partition function of spinless CS model from that of the spin CS 
model. By using this partition function of SU (m) HS spin chain, it is possible to study 
the energy level density distribution and the nearest neighbour spacing distribution 



such HS spin- 2 chain with long-range interaction was originally motivated from the 




(1.1) 
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for fairly large values of N [9]. Interestingly, it has been found that, the continuous 
part of such energy level density follows Gaussian distribution to a high degree of 
approximation. 

In this context it may be noted that, there exists a SU{'m\n) supersymmetric ex- 
tension of HS spin chain [3], where each site is occupied by either one of the m type 
of bosonic states or one of the n type of fermionic states. Such supersymmetric spin 
chains play an important role in describing some correlated systems of condensed mat- 
ter physics, where holes moving in the dynamical background of spin behave as bosons 
and spin- 1 electrons behave as fermions [13]. It is worth noting that the supersym- 
metric 5'C/(m|n) HS spin chain exhibits Y{gl(^m\n)) super- Yangian symmetry [3], which 
is also the quantum group symmetry of supersymmetric SU{m\n) Polychronakos spin 
chain [14,15]. Consequently, by using the motif representations and skew- Young di- 
agrammes associated with supersymmetric Polychronakos spin chain [15], one can in 
principle calculate the degeneracy factors for all energy eigenvalues of supersymmet- 
ric HS spin chain. However, similar to the nonsupersymmetric case, this method for 
finding the full spectrum and related partition function becomes very complicated for 
large values of N. 

The aim of the present article is to find out the exact partition function for su- 
persymmetric SU{m.\n) HS model by applying the freezing trick and also to study 
global properties like level density distribution of the corresponding spectrum. For this 
purpose, it is convenient to map the supersymmetric HS model to a usual spin chain 
containing an 'anyon like' representation of the permutation algebra as spin dependent 
interactions [16,17]. In Sec. 2 we describe this mapping and also show how the freezing 
trick can be applied for the case of 5'C/(m|n) HS spin chain by embedding it in a spin 
CS model containing the same anyon like representation of the permutation algebra. 
In Sec. 3, we find out the complete spectrum of such spin CS model including the de- 
generacy factors for all energy levels. In Sec.4, we calculate the partition function of 
this spin CS model at the strong coupling limit and divide it by that of the spinless 
CS model to finally obtain the partition function of SU{m\n) HS spin chain. In this 
section, we also discuss about the motif representation for SU{m\n) HS spin chain 
and find that, due to the lifting of a selection rule, some extra energy levels appear 
in the spectrum in comparison with the case of SU (m) spin chain. Subsequently, we 
conjecture that the partition function of SU (m|n) HS model satisfies a duality relation 
under the exchange of bosonic and fermionic spin degrees of freedom. In Sec.5, we 
study the level density distribution and the nearest neighbour spacing distribution for 
the spectrum of 5'C/(m|n) HS spin chain by using its exact partition function. It is 
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found that, for sufficiently large values of N, continuous part of the energy level den- 
sity satisfies the Gaussian distribution with a high degree of accuracy. We also derive 
exact expressions for the mean value and standard deviation which characterize such 
Gaussian distribution. Sec. 6 is the concluding section. 

2 Application of the freezing trick 

For the purpose of defining the SU (m|n) supersymmetric HS spin chain, let us consider 
a set of operators like Cj^{Cja) which creates (annihilates) a particle of species a on 
the j-th lattice site. These creation (annihilation) operators are assumed to be bosonic 
when a e [1, 2, m] and fermionic when a e [m + 1, m + 2, m + n]. Thus, the 
parity of C't (C'ja) is defined as 



p(Q«) = p(Cjj = Oforae [1,2, ....,m] , 

p{Cja) = p{Cja) = 1 for a e [m + 1, m + 2, m + 



These operators satisfy commutation (anti-commutation) relations hke 

[Cja, Ckp]± = , [Cj^, Clf^]± = , [Cja, Clf^]± = SjkSap , (2.1) 

where [A, B]± = AB - (-l)P(^)p(^)fiA. Next, we focus our attention to a subspace of 
the related Fock space, for which the total number of particles per site is always one: 

m+n 

J2 ClCja = 1, (2.2) 

a=l 

for all j. On the above mentioned subspace, one can define supersymmetric exchange 
operators as 

m+n 
a,/3=l 

where 1 < j < k < N. These -Pjj^'"^'s yield a realization of the permutation algebra 
given by 

V%^1, VjkVki = VjiVjk = VkiVji , [Pjk, Vim] = , (2.4) 

where j. A;, Z, m are all distinct indices. Replacing Pjk by -Pj^'"^ in eqn.(l.l), one 
obtains the Hamiltonian of SU{m\n) supersymmetric HS model as [3] 



( 1 + P "'"M 



2 ^ sin^f^,- - ik) 

l<j<k<N ^'^^ 



4 



Now we want to describe how this 5'C/(m|n) supersymmetric HS model (2.5), con- 
taining creation-annihilation operators, can be transformed to a spin chain. To this 
end, we consider a particular type of anyon like representation of permutation algebra 
(2.4), which acts on a spin state like \a1a2 ■ ■ ■ ajv) (with aj e [1,2,..., m-|-n]) as [16,17] 

Pj™l"Vi . . . a,- . . . afe . . . ajv) = e'*(">'"^+i'-'"'=)|ai ■ ■ ■ . . . a,- . . . ajv), (2.6) 

where e^*("^'"^+i'-'°'=) = 1 if a^, 0;^ G [1, 2, . . . , m], e'^^^^'^^+i--'"*) = -1 if a^-, e [m-F 
l,m + 2, . . .,m + n], and g^^W'^^+i'-'^*) = (-l)'^^P=i+i^r=m+i*"p.- jf e [1,2, ... ,m] 
and ctfc e [m + 1,171 + 2, . . . ,m + n] or vice versa. For the purpose of interpreting the 
phase factor e»*{o=3."j+i. ^ physical way, it is convenient to call a, a 'bosonic' spin 
when ai G [1,2,..., m] and a 'fermionic' spin when e [m -|- 1, m + 2, . . . , m -|- n]. 
Prom eqn.(2.6) it follows that, the exchange of two bosonic (fermionic) spins produces 
a phase factor of 1 (—1) irrespective of the nature of spins situated in between the j-th 
and fc-th lattice sites. However, if we exchange one bosonic spin with one fermionic 
spin, then the phase factor becomes (—1)'^ where p is the total number of fermionic 
spins situated in between the j-th and k-ih lattice sites. 

Next we observe that, due to the constraint (2.2), the Hilbert space associated 
with SU{m\n) HS Hamiltonian (2.5) can be spanned through the following orthonor- 
mal basis vectors: Cl^_^Cl^^ . . . Clj^^\0) , where |0) is the vacuum state and aj e 
[1, 2, . . . , m + n]. Consequently, it is possible to define a one-to-one mapping between 
these basis vectors and those of the above mentioned spin chain as 

\a,a2 ...aj,)^ ClA. ■ ■ ■ ^kjO). (2.7) 

With the help of commutation (anti-commutation) relations (2.1), one can easily verify 
that 

p(m|n) /^t r't r't in\ 

= e^*(".v-'"'=) ...Cl^... Cl^ . . . CLJO), (2.8) 

where e**^"^' "'"*^ is the same phase factor which appeared in eqn.(2.6). Comparison 
of eqn.(2.8) with eqn.(2.6) through the mapping (2.7) reveals that, the anyon hke 
representation Pj^'"^ is equivalent to the supersymmetric exchange operator Pj^'"^ 
Hence, if we define a spin chain Hamiltonian through Pj^'"^ as 

p{ni\n) 



TT{m\n) _ 1 V ' (2 Q) 



2 ^ sin^(^. - ik) 

l<j<k<N ^'^^ 
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that would be completely equivalent to the supersymmetric SU{m\n) HS model (2.5) 
[16]. Clearly, for the special case n — 0, -Pjj^'"^ reproduces the original spin exchange 
operator Pj^ and (2.9) reduces to the Hamiltonian of SU (m) HS spin chain (1.1). 

Since it is convenient to apply the freezing trick to the spin chain Hamiltonian (2.9), 
for the rest of this article we shall deal with this form of supersymmetric SU{m\ri) HS 
model instead of its original form (2.5). 

By using the anyon like representation -Pjj^'"\ one can construct a spin CS model 

like 

which contains spin as well as dynamical degrees of freedom and the positive parameter 
a as coupling constant. With the help of mapping (2.7) it can be shown that, this spin 
CS model is equivalent to a supersymmetric spin CS model [18] with Y{gl(^^\n)) super- 
Yangian symmetry. The spin CS Hamiltonian H* (2.10) might be formally written 
as 

//* = //o + 4aH("*l"), (2.11) 
where Hq is the Hamiltonian of spinless CS model given by [19] 

j=l ^ ^<J<k<N ^ ■' 

and H*^™l"^ is obtained from if^*"!"-) (2.9) by the replacement Xj. Now the decou- 

pling of the dynamical degrees of freedom of H* (2.10) from its spin degrees of freedom 
can be achieved by using the freezing trick [10-12]. This trick is based on the fact 
that in the Umit a — > oo, particles freeze at the equilibrium positions of Hq, which are 
simply the lattice points {^j) of the spin chain in eqn.(2.9). Consequently, by using 
eqn.(2.11) at freezing limit, we find that the energy levels of H* are approximately 
given by 

E*^^Eoj + AaEk, (2.13) 

where Eqj and Ek are any two levels of Hq and respectively. Hence, we obtain 

a relation like 

Z(-I")(T) = lim (2.14) 

where Z^"^^"'\ Z* and Zq denote the partition functions corresponding to the Hamil- 
tonians H* and Hq respectively. Thus the freezing trick allows us to compute 

the partition function of SU (m|n) supersymmetric HS spin chain, by modding out the 
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contribution of spinless CS model from the partition function of spin CS model (2.10). 
Due to the Gallehan invariance of H* and Hq it follows that, if ^ is an eigenstate of 
any one of these Hamiltonians with momentum p, then ■0' = e^''^^3=i^j'0 will also be 
an eigenstate of the same Hamiltonian with momentum {p-\-2tN). As a result, we can 
always adjust the parameter r such that ■0' will be an eigenfunction of H* or Hq with 
zero momentum. In this article, we shall always consider eigenstates of these Hamilto- 
nians with zero momenta and evaluate the partition functions Z* as well as Zq at the 
center of mass frame. Since both Z* and Zq get modified by the same multiplicative 
factor due to a Gallelian transformation, does not depend on the choice of the 

reference frame. 



3 Spectrum of spin CS model 

In this section our aim is to find out the complete spectrum of spin CS model (2.10) 
containing anyon like representation of the permutation algebra. Even though the 
spectrum of such spin CS model has been studied earlier [17], multiplicities of degen- 
erate eigenfunctions corresponding to all energy levels have not been found. Since 
these numbers are required for calculating the partition function of this model, here 
we want to derive a general expression for the degeneracy factors of all energy levels. 
It is well known that the eigenfunctions of spin CS Hamiltonian (2.10) can be written 
in a factorised form like 

il){xi, . . . , xat; ai, . . . , Oat) = r''0(a;i, . . . , xn] ai, . . . , q;2), (3.1) 

where F = ni<7 sin(a;i — Xj). By operating H* (2.10) on the above form of ■0, we find 

= r"^>, (3.2) 



that 



where 



H* ^4 



d \ ^ sr-^ Zk + Zj f d d 



Zk r\ Zj - 



dZj J k<' — ^ 

k<j \ J / 



(3.3) 



with Zj = e^*^^ . Equation (3.2) implies that, if is an eigenvector of H* with eigenvalue 
E, then would be an eigenvector of H* with the same eigenvalue. Thus the 
diagonalisation problem of H* is reduced to the diagonalisation problem of H*. 
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For solving H*, it is convenient to introduce another operator H which acts only 
on the coordinate degree of freedom and may be given by [7] 



n 



k<j \ 3 Kj 



(3.4) 



where the Kj^ is the coordinate exchange operator which exchanges the coordinates of 
j-th and /c-th particle: 

Kjkf{xi, . . . , Xj, . . . , Xk, . . . , Xjv) — f{xi, . . . , Xk, . . . , Xj, . . . , xn). 

It may be observed that, H* (3.3) can be reproduced from the expression of Ti (3.4) 
through the substitution Kjk — Pj™'"^ This connection between H* and Ti will play 
a crucial role in our calculation for finding the spectrum of H*. Let us first consider 
state vectors given by monomials like 

— z\^z^ ■ ■ ■ z^N ^ (3-5) 

where p = {pi,P2, ■ ■ ■ ,Pn} £ satisfies the constraints: (i) {pi — Pj) are integers 
for all i,j and (ii) Y^f^iPi — 0- The last condition implies that these monomials 
represent state vectors with zero total momentum. In particular, one can consider 
corresponding to a nonincreasing vector p = (pi,p2, . . . ,Pn), whose elements satisfy 
the conditions: (i) li = pi — pi+i is a nonnegative integer for i e [1, . . . , A^" — 1] and 
(ii) J2iLi Pi — 0- is evident that, N — 1 number of nonnegative integers (Zj's) are 
sufficient to specify a nonincreasing vector p. Given two distinct nonincreasing vectors 
p and p', we shall write p ^ p' if pi — p'^ = • • • = —Pi_i — and pi < p\. A partial 
ordering can be defined on monomials like (3.5) in the following way. By permuting 
the elements of any p, one can always construct a unique nonincreasing vector p. The 
basis element would precede ^p> if p ^ p', where p and p' are nonincreasing vectors 
obtained from p and p' respectively by permuting their components. The above defined 
ordering is effectively a partial ordering, since it does not induce an ordering between 
^p and ^p/ when p = p'. It can be shown that the action of H (3.4) on the state vector 
Cp yields [7,9] 

mp = Ep^p+ Y] cpp<p', (3.6) 



p' 
(p'<p) 



where 

N 



Ep = Y,{'2Pi + a{N + l-2i)}\ (3.7) 



i=l 
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Thus it is clear that, if one constructs a Hilbert space through basis vectors of the 
form (3.5) and partially order them in the above mentioned way, then Ti. will act as a 
triangular matrix on this space. 

Next we want to construct another partially ordered Hilbert space, on which H* 
(3.3) can be represented as a triangular matrix. To this end, we define a set of permu- 
tation operators as n^-^'"-' = P^'^^'^^Kj^. Since both -Pj^'"^ and Kjf. satisfy an algebra of 
the form (2.4), while acting on the spin and coordinate spaces respectively, the newly 
defined operator H^™'"^ also yields a representation of the same permutation algebra 
on the direct product of coordinate and spin spaces. Hence, by using this representa- 
tion of permutation algebra, we can construct a 'generalized' antisymmetric projection 
operator A^"*'") satisfying the relation 

]^KI")^(m|n) ^ ^{m\n)jj{Tn\n) ^ _^(m|n) /g^g-v 
jk jk ' V / 

or, equivalently, pj^l")A("»l") = -XjfeA("»l") [16,17]. Even though A("^l") can be ex- 
pressed as a function of H^-^'"\ explicit form of this projection operator is not neces- 
sary for our present purpose. However it may be noted that, since both Kjk and -Pj^'"^ 
commute with H (3.4), A*^"^l") also satisfies the relation 

[n, A('"l'^)] = . (3.9) 

With the help of projection operator A("*l"\ we define a state vector on the direct 
product of coordinate and spin spaces as 

0p = " = A^"^'"Hepl«i ■ ■ ■ ocn)}. (3.10) 
Using eqns.(3.8) and (2.6) it can be shown that 

<Kl::Z:Z:Z = -A^"^'"^i^.fc^r'"^{^f ... .. ... 4^ |ai a^)} 

= _e^*K v."^)A(-l"){zf . . . 4* . . . ^-^ . . . 4- |ai . . . ftfc . . . a, . . . a^)} 

— _„i^{aj,...,ak)±ai...ak---aj...aN 

- e <Ppi...p^...pj...pt, ■ l<3.iij 
By repeatedly using the above equation we find that 

0^ = e(a,p)0?', (3.12) 

where e{ot, p) = ±1, p is the nonincreasing vector corresponding to p and oc' is a spin 
vector which is obtained by permuting the components of a. Hence, all state vectors 
of the form (3.10) can be obtained by choosing p from the set of nonincreasing vectors 
only. 
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Corresponding to any given nonincreasing vector p, one can define a vector space 

as 

^p={<f>p \ai,...,aN e [1,2,... m + n]). (3.13) 

It is important to note tfiat, different values of a. may lead to the same 0? which is 
a basis element of V^. To see this thing in a simple way, let us take a nonincreasing 
sequence p satisfying the condition pi = Pj = p (say). For this special case, eqn.(3.11) 
reduces to 

Vp^...p...p...p^ c rpi...p...p...p,^ ■ l-J'-L^; 

Clearly (l>^\ ""^p "^^' pl^ and 4>'p\ '"p "'p ' 'pl^ represent the same state vector (up to a phase 
factor), although they correspond to different values of a. For a given (j)p, we say 
that two spin components of o: belong to the same 'sector' if the corresponding two 
components of p are equal to each other. For example, the spin components ai and 
a,- appearing in the state 0"! '"^"'"^ belong to the same sector according to this 
convention. Since e** = 1, for ai,aj G [1,2, ...,m], it is clear from eqn.(3.14) that 
bosonic spins within the same sector obey 'fcrmionic statistics' after antisymmctrisa- 
tion. In particular, two bosonic spins of same flavour can not coexist within a single 
sector. Similarly, since e** = —1 for Oj, G [m + 1, + 2, . . . , m + n], one can flnd 
from cqn.(3.14) that fcrmionic spins within the same sector obey 'bosonic statistics' 
after antisymmctrisation. Therefore, any number of fcrmionic spins having the same 
flavour can be accommodated within a single sector. 

Now we want to flnd out the dimensionality of the space V^. For this purpose, it 
is useful to write p in the form 

ki ki kr 

p = (pi7~^^Tpi, • • • ,pir''^^TA, • • • (3-15) 

where pi > • • • > pi > ■ ■ ■ > Pr, Yll=i ^« — integer which can take 

any value from 1 to N. It is obvious that k = {k-y, . . . ,kr}, which belongs to the 
set Vn of ordered partitions of N, may be treated as a function of p. For a given 
(f)*^, clearly the components of a are separated into r different sectors where the i-th 
sector contains ki number of spins. It is evident that the dimensionality of the space 
Yp may be obtained by counting the number of independent ways one can distribute 
total N number of spins within r sectors. To this end, let us first try to find out the 
number of independent ways of filling up the i-th sector through ji number of bosonic 
spins and j2 number of fermionic spins, where ji + j2 = k^. Using eqn.(3.14) we have 
already seen that, bosonic and fermionic spins within the same sector obey fermionic 
and bosonic statistics respectively. Therefore, we can pick up ji number of bosonic 
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spins from m different flavours in '^Cj^ different ways and j2 number of fermionic spins 
from n different flavours in ^'^'^'^'^Cj^ ways, where ^Ci — ^, for Z < p and ^Ci — 
for I > p. Thus the number of independent ways of filling up the i-th sector through 
ji number of bosonic spins and j2 number of fermionic spins is given by 



Summing up these numbers for all possible values of ji and j2, we obtain the total 
number of independent ways of filling up the i-th sector through ki number of spins as 

ji= 

min (m,fei) , . , 

V- m! {n + kj-ji-iy. 



Since two spins belonging to different sectors do not follow any exchange relation like 
(3.14), the number of independent ways we can distribute total number of spins 
within r different sectors is given by the product of all d^^^'^^{ki). Therefore, by using 
(3.16), we finally obtain the dimension of as 

dH")(k) = JJ(i("^l")(A;0 = n "^C,-"+'=*-^-'Cfc._J . (3.17) 

i=i i=i y j=o j 

Even though this expression is derived by assuming that bosonic and fermionic spin 
degrees of freedom (i.e., m and n respectively) take nonzero values, it is also possible 
to obtain the dimension of Vp for the SU{n) fermionic case by putting m = in 
eqn.(3.17): 

r 

c^(0N(k) = -+*=>-iCfc,. (3.18) 

i=l 

Furthermore, by putting n = in eqn.(3.17), subsequently using the relation = 
for Z > p > 0, and also assuming that "^Cq = 1, one can reproduce the dimension of 
Vp for the SU (m) bosonic case [9] as 

r 

dH0)(k) = -Q mc^^_ (3.19) 

i=l 

It is interesting to observe that, while d*^"^l°)(k) (3.19) can take a nonzero value only 
if A;^ < m for all i, both d("^l")(k) (3.17) and d(o|")(k) (3.18) take nonzero values for 
any k e Vn- Consequently, Vp will represent a nontrivial vector space for the SU{m) 
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bosonic case only if at most m components of p take the same value. On the other 
hand, Vp will represent a nontrivial vector space for all possible values of p when at 
least one fermionic spin degrees of freedom is present. 

The Hilbert space associated with H* (3.3) may now be defined by taking the direct 
sum of Vp (3.13) for all allowed values of p : 

V = ©V^. (3.20) 
p 

We define a partial ordering in this Hilbert space by saying that the basis element (pp 
precedes 0^' if p -< p'. By consecutively applying the relations (3.10), (3.8), (3.9), 
(3.6) and (3.12), it is easy to check that 

^*(0?) = AH-)7i^.|Q,i...a^) 

= a("^1^)(ep^p+ J2 Cpp'ep')|ai---aiv) 

^ P' ^ 

(P'<P) 

= Ep<j>^+ J] e(a,p')c^p'0^;. (3.21) 

P' 
(p'<p) 

Hence H* is represented as a triangular matrix on V. Diagonal elements of this trian- 
gular matrix yield the eigenvalues of H* as 

N 

E*{p, cx) = Ep^Yl + a{N + l- 2i)}'. (3.22) 

i=l 

Consequently, the eigenvalues of spin CS Hamiltonian H* (2.10) are also given by 
E*{p, a) in the above equation. 

Since E*{p, ex.) in eqn.(3.22) does not really depend on the spin vector a, the number 
of degenerate energy eigenstates associated with the quantum number p would coincide 
with the dimension of the space V^. Thus the degeneracy factor of the energy eigenvalue 
E*{p,cx.) corresponding to the quantum number p is given by (i("*l")(k) appearmg m 
eqn.(3.17). We have already seen that, in contrast to the pure bosonic case, d^"*l")(k) 
takes nonzero values for all possible k e Vn when at least one fermionic spin degrees of 
freedom is present. Consequently, the presence of fermionic spin degrees of freedom in 
H* (2.10) would lead to a spectrum with many additional energy levels in comparison 
with the spectrum of bosonic spin CS model. 

Finally let us briefly comment about the known spectrum of spinless CS Hamilto- 
nian Ho (2.12) [19]. Using the fact that the eigenfunctions of Hq can be written in a 
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factorised form like ipo — ^"'(poi^i, ■ ■ ■ , ^n), it is possible to transform Hq into Hq as 

i/o^o = ™o0o, 

where can be obtained from H (3.4) through the substitution Kij — > 1. For 
constructing the Hilbert space associated with one may consider elements like 
0p = Ao(^p), where Aq is the symmetriser in the coordinate space: Kj^Aq — Aq. Since 
0p = where p is the nonincreasing vector corresponding to p, the Hilbert space of 
^0 is defined through independent basis vectors (j)p for all values of p. An ordering can 
be defined among these state vectors by saying that 0p precedes 0p/ if p ^ p'. Using 
eqn.(3.6) it can be shown that, Hq acts as a triangular matrix on these completely 
ordered basis vectors and the eigenvalues of Hq are also given by Ep in eqn.(3.22). 
However, due to the absence of spin degrees of freedom, only one energy eigenstate is 
obtained corresponding to each quantum number p in this case. 

4 Partition function of SU{m\n) HS spin chain 

By using the freezing trick we have seen that, the partition function of supersymmetric 
SU{m\n) HS spin chain can be obtained by dividing the partition function of spin 
CS model (2.10) at the strong coupling limit through that of the spinless CS model 
(2.12). To execute this programme, let us first briefly recapitulate the calculation for 
the partition function of spinless CS model (2.12) at a ^ oo limit [9]. It should be 
noted that, the eigenvalues in eqn.(3.22) can be expanded in powers of a as 

TV 

E*{p, a) = E{p) = a^Eo + 4a ^(iV + 1 - 2t)pi + 0(1), (4.1) 

i=l 

where Eq = ^N{N'^ — 1) . Since Eq does not depend on p or a, the effect of this Eq 
will be manifested as the same overall multiplicative factor in the partition functions 
of spin CS model and its spinless counterpart. Hence, by dropping the first term in 
eqn.(4.1), and neglecting the 0(1) term in the limit a — > oo, one can write down the 
partition function of spinless CS model (2.12) as 

Zo(4ar) ~ q^'^'^''^'-^'^ , (4.2) 
p 

where q — e~^^^'"^^\ Using N — 1 number of nonnegative integers (/j's) which uniquely 
determine p, one can further simplify this partition function as [9] 

N-l N-1 ^ 

z„(4aT)c. y: n n (^-s) 

h,...,lN_l>0 j=l j=l 
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Next, we want to calculate the partition function of spin CS Hamiltonian (2.10) at 
a — > oo limit. Dropping again the first term as well as 0{1) term from the right hand 
side of expansion (4.1), and expressing the nonincreasing vector p through eqn.(3.15), 
£'*(p, a) can be written as 

r Ki 

£;*(p,a)~4a^p, ^ {N + l-2j), (4.4) 

i=l j=Ki_i+l 

where Ki = Yl]=i denote the partial sums corresponding to the partition k G Vn 
and Kq = 0. Using a set of variables like Uj = pj — pj+i for j G [1, . . . ,r — 1] (since 
Pj > Pj+i, all tA,'s are positive integers), one can express the energy eigenvalue in 
eqn.(4.4) as 

r-l 

E*{p,cx)c^4aJ2^jNj, (4.5) 

i=i 

where Nj = Kj{N — Kj). It may be noted that, due to the condition X^i^iP* ~ 0; 
r — 1 number of z/j's uniquely determine the nonincreasing vector p in eqn.(3.15). Con- 
sequently, the single sum can be replaced by the double sum Vn '^^i ur-i>o 
in the expression of the partition function. By using the eigenvalue relation (4.5) and 
the corresponding degeneracy factor ci('"l")(k) (3.17), we obtain the partition function 
of spin CS Hamiltonian (2.10) at a — > oo limit as 

r-l 

Z*{AaT) ~ ^(i('"l")(k) JJg^^'^^- 

ke Vn vi,...,Vr-i>Q 3=1 

r-l 



r-l TV- 

E ^^'"'"^(k)ni4^- (4-6) 



Q 

ke Vn J=1 ^ 



Using eqns.(2.14), (4.3) and (4.6), we finally obtain the partition function of 5'[/(m|n) 
HS spin chain as 



Af-l r-l 



qN, 



1=1 ke Vn i=i ^ ^ 

Since the partial sums Ki, K2, . . . , Kr-i associated with k are natural numbers obeying 
1 < Ki < ■ ■ ■ < < N — 1, one can define their complements {K'- 's) as elements of 
the set: {1, . . . , N — 1} — {Ki, . . . , Kj.-i}, which satisfy the ordering K[ < • ■ ■ < K'j^_^. 
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Hence one can rearrange the product nili^(l ~ i^^^ terms as [9] 

JV-l r-1 N-r 

11(1-?'^"-'^)= 11(1 -^"011(1-^"')' (4.8) 

1=1 j=l i=l 

where — K-{N — K'^. By substituting this relation to eqn.(4.7), we get a simphfied 
expression for the partition function of SU{m\n) HS model as 

^Hn)(y) _ ^ ^(m|n)(j^) ^,=1 Y[{l - q""'^) . (4.9) 
ke Vn «=1 

We have already seen that, both d("^l")(k) (3.17) and d(°l")(k) (3.18) take nonzero values 
for any k e V^sf- Consequently, in contrast to the restricted choice of k for the case 
of bosonic spin chain [9] , all possible k e will contribute to the partition function 
(4.9) in the case of supersymmetric as well as fermionic HS spin chain. 

It is well known that the spectrum of bosonic SU (m) HS spin chain (1.1) containing 
N number of lattice sites can be obtained from motifs like 5 = (0, 5i, . . . , 5n-i-i 0), where 
each 5j is either or 1 [6-8]. The form of these motifs and corresponding eigenvalues 
can be reproduced by using the partition function of bosonic SU{m) HS spin chain 
[9]. Now we want to explore how the motifs associated with SU{m\n) supersymmetric 
HS spin chain emerge naturally from the expression of partition function (4.9). To 
this end, we define a motif corresponding to the partition k by using the following 
rule: 5j = if j coincides with one of the partial sums Ki and 5j — 1 otherwise. 
Furthermore, it is assumed that the lowest power of q in eqn.(4.9) for the partition k 
gives the energy eigenvalue E{S) of the above motif S. In this way we obtain the energy 
levels of SU{m\n) supersymmetric HS spin chain as 

r—l J.J-, j.j-2 \ N—1 

E{S) = iV, = ^ ~ ^ + J2 - N) , (4.10) 

i=i j=i 

which apparently coincides with that of the bosonic HS spin chain. However it should 
be observed that, for the case of SU{m) spin chain, only those k would contribute in 
the partition function for which Kj — Kj_i — kj < m [9]. This leads to a selection rule 
which prohibits the occurrence of m or more consecutive I's within the corresponding 
motifs. On the other hand, since all k e Pn contribute to the partition function (4.9) 
of supersymmetric HS spin chain, it is possible to place any number of consecutive I's 
or O's within a motif 5. Consequently, the selection rule occurring in the bosonic case 
is lifted for the case of supersymmetric HS spin chain and many extra energy levels 
appear in the corresponding spectrum. This absence of selection rule in the spectrum 
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of supersymmetric HS spin chain was previously observed by Haldane on the basis 
of numerical calculations [3]. By using the expression of E{5) in eqn.(4.10), we can 
easily evaluate the maximum and minimum energy eigenvalues of this system. Prom 
the expression of E{S) it is evident that, the motif S = (0, 0, . . . , 0, 0) would correspond 
to the maximum energy Emax — '^^'^q ■ Similarly for the motif 5 = (0, 1, . . . , 1, 0), 
we obtain the minimum energy of the system as Emin — + X^jlJi^ i(i — N) — 

0. It is interesting to note that these maximum and minimum energy eigenvalues of 
SU{m\n) supersymmetric HS spin chain do not depend on the values of m and n. 
Moreover, the lifting of the selection rule is responsible for the zero minimum energy 
of supersymmetric HS spin chain. 

Using Mathematica we find that, for a wide range of values of m, n and N, the 
partition function (4.9) of SU{m\ri) HS model satisfies a duality relation of the form 

Z("*l")(g) = g^^Z("l'")(g-^). (4.11) 

This result motivates us to conjecture that the above duality relation, involving the 
interchange of bosonic and fermionic spin degrees of freedom, is valid for all possi- 
ble values of m, n and N. It may be noted that, for the particular case n — 0, 
eqn.(4.11) relates the partition function of SU{m) bosonic HS spin chain to that of 
SU (m) fermionic spin chain. By applying the relation Pj^^^^ — —Pj^^^ and the sum- 
mation formula ^i<j<fc<jv sin^(g'^._g^) = ^^^^ [9,20], we find that the Hamiltonians 
of bosonic and fermionic spin chains are connected as 

^(m|0) ^ ^(^^ ~ ^) _ ^(0|m) 

6 

Using the above relation along with the definition of partition function given by 

one can easily prove eqn.(4.11) for the particular case n = 0. It 
would be interesting to explore whether eqn.(4.11) can be also proved for the general 
case by establishing some relation between Pj™'"'' and -Pj^'™^ Comparing the coeffi- 
cients of same power of q from both sides of eqn.(4.11), we find that the energy levels 
of SU{n\m) spin chain can be obtained from those of SU{m\n) spin chain through the 
transformation Ei — > ^'^^q — Ei and also get the relation 

p(m|n) ^^.^ ^ p(n|m) ^ ^(^^ ~ ^) _ E^j , (4.12) 

where V^'^\'^\Ei) denotes the degeneracy factor corresponding to energy Ei of SU {m\n) 
HS spin chain. Thus it is evident that, the spectrum of SU{n\m) spin chain can be 
obtained from that of SU{m\n) spin chain through an inversion and overall shift of all 
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energy levels. Such relation between the spectra of supersymmetric HS spin chains was 
empirically found by Haldane with the help of numerical analysis [3]. 

5 Spectral properties of SU{m\n) HS spin chain 

In this section we shall explore some spectral properties of supersymmetric HS model 
by using its exact partition function ^("^'"■^(T') (4.9). It has been already mentioned 
that, calculation for the degeneracy factors associated with the energy eigenvalues of 
this spin chain becomes very cumbersome by using the motif representations for large 
values of N. However, with the help of a symbolic software package hke Mathematica, 
it is possible to express the partition function (4.9) as a polynomial of q and explicitly 
find out the degeneracy factors of all energy levels for relatively large values oi N. In 
this way, we can study properties like level density distribution and nearest-neighbour 
spacing (NNS) distribution for the spectrum of supersymmetric HS spin chain. 

For the case of SU (m) bosonic spin chain, it has been found earlier that the con- 
tinuous part of the energy level density obeys Gaussian distribution to a very high 
degree of accuracy ior N » 1 [9]. At present, our aim is to study the level density 
distribution in the spectrum of SU (m|n) supersymmetric HS spin chain and investigate 
whether it exhibits a similar behaviour. To begin with, let us consider the simplest 
case of SU supersymmetric HS spin chain. In this case, the degeneracy factor 
in eqn.(3.17) reduces to a simple form given by (i(^l^)(k) = 2^ By substituting this 
degeneracy factor to eqn.(4.9), taking some specific value for the number of lattice 
sites hke N — 15 and using Mathematica, we express the partition function of SU 
spin chain as a polynomial of q. The coefficient of in such polynomial evidently 
gives the degeneracy factor X'(^l^)(£'j) corresponding to the energy eigenvalue Ei, which 
we plot in Fig.l. This figure clearly indicates that the energy level distribution obey 
Gaussian approximation but with some local fiuctuations. Similar behaviour of energy 
level distribution has been found by studying SU (m|n) HS spin chain with other values 
of m, n and sufficiently large values of N. 

Prom the above discussion it is apparent that, if we decompose the energy level 
density associated with SU{'m\n) HS spin chain as a sum of continuous part and fiuc- 
tuating part, the continuous part will obey Gaussian distribution for large values of 
N. This behaviour of the continuous part can be measured in a quantitative way by 
studying the cumulative level density [9], which eliminates the fluctuating part of the 
level density distribution. For the case of SU{'m\n) HS spin chain, cumulative level 
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density of the spectrum is defined as 



F{E) 



(m + n) 



N 



(5.1) 



Ei<E 



Obviously, this F{E) can also be obtained by expressing the exact partition function 
(4.9) as a polynomial of q. We want to check whether this F[E) agrees well with the 
error function given by 



1 + erf 



E- n 



(5.2) 



where and a are respectively the mean value and the standard deviation associ- 
ated with the energy level density distribution. These parameters are related to the 
Hamiltonian ii"^'"!'*) (2.9) as 



tr [//("^l")] 



wv ' = , ^ .N /^'- (5.3a, 6) 

[m + n) [m + n) 

For the purpose of comparing F{E) with G{E), it is necessary to express the pa- 
rameters // and (T as some functions of m, n and N . To this end, we need the following 
trace formulas: 



tr 
tr 
tr 



{m\n)\2 
ij ) 



tr [1] 



tr 



P. 



{m\n) 



f){m\n) f){m\n) 

Y){m\n) f3(™l") 
^ij ^kl 



tr 

N-A,2 



f){m\n) p(m|n) 



tr 



p{m\n) p(m\n) 
ij kj 



(5.4a, 6) 
(5.4c) 

(5.4rf) 



where s = m + n, t — m — n and k, I are all different indices. Derivation of these 
trace formulas is given in Appendix A of this article. For obtaining the functional form 
of II and cr, it is also required to evaluate summations like 



-^0 = V] . 2/e ' P^ = ■ 4/t T^- 

^ sm - Q ^ sm - Cj) 



i<j k<l 



sin^(Ci -Cj)sin^(Cfe -6) ' 



(5.5a, 6) 

(5.5c) 



(5.5ci) 



It is easy to see that the above defined Rq, Ri, R2 and R3 satisfy the relation 

Rl^Ri + R2 + R3. (5.6) 
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Using some summation formulas given in Ref.20, it can be shown that 

R.^!^^, (5.7a) 

(5.7.) 

_ 7V(iV^ - 1)"(7V - 4) 7^(7^^= - l)(jV^ + 11) 
^ 36 90 ' ^ ' 

^^^^Nl^ml^ (5.7d) 

Derivation of these relations is discussed in Appendix B of this article. 

Now, by using eqns.(5.3a), (5.4a,b) and (5.7a), we can express /i. as a function of 
m, n and N given by 

Next, by using the trace formulas (5.4a,b,c,d), we obtain 

oN~2( 2,nf\ N ^N-4f2 N-2 

tr [{Hi-\-^r] = ^ Rl+'-^R. + ^ R2+'-^Rs. (5.9) 

Substituting the expressions for fi in eqn. (5.8) andtr [(i7(Hn))2j 

in eqn. (5.9) to eqn. (5.3b), 

and subsequently using (5.6), it can be shown that 

s4 _ /2 2 _ ^2 

Finally, by substituting the values of Ri (5.7b) and Rs (5.7d) to eqn. (5. 10), we can 
express cr as a function of m, n and N given by 

'—±N(N'-l)(N' + n) + '-^N(N'-l){N'-4)y . (5.11) 

Thus we are able to find out the functional forms of the parameters /i and a for the 
case of SU{m\ri) HS spin chain. It may be observed that, in the special case n — 
(for which one gets s — t — m), eqns.(5.8) and (5.11) reproduce the forms of /i and a 
corresponding to the SU (m) bosonic HS spin chain [9] . 

Now we can compare the cumulative level density F{E) (5.1) with the error func- 
tion G{E) (5.2), where the values of // and a are obtained from eqns.(5.8) and (5.11) 
respectively for any given m, n and N. In Fig. 2, we plot such F(E) and G{E) for 
the particular case of SU spin chain with A^" = 15 lattice sites. From this figure 
it is evident that F{E) follows G{E) to a high degree of approximation. One can 
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also quantify the agreement between F{E) and G{E) by calculating the corresponding 
mean square error (MSE), which for the above mentioned case is given by 8.46 x 10~^. 
It may be noted that, the agreement between F{E) and G{E) improves rapidly with 
increasing values of N. For example, in the case of SU{1\1) model, MSE between 
F{E) and G{E) decreases from 5.17 x 10"^ to 8.46 x 10"^ when the value of N is 
increased from 10 to 15. Next, we consider the particular cases of SU {1\2) as well 
as SU {2\1) supersymmetric spin chain with N — 15 lattice sites and also the SU{3) 
bosonic spin chain with same number of lattice sites for the sake of comparison. In 
Fig.3, we plot F(E) and G(E) for SU (1\2), SU (2\1) and SU(3) HS spin chains and 
find that the corresponding MSEs are given by 2.38 x 10-^ 1.2 x 10~^ and 8.61 x 10"^ 
respectively. Again F{E) shows very good agreement with G{E) for all of these cases. 
Such agreement also improves rapidly with increasing values of N. For example, in the 
case of SU{1\2) spin chain, MSE between F{E) and G{E) decreases from 8.23 x 10"^ 
to 2.38 X 10~^ when the value of N is increased from 10 to 15. Analysing many other 
particular cases with different values oim, n and sufficiently large values oi N, we find 
that F{E) follows G{E) with a high degree of approximation for all of these cases. 

From the above discussion it is evident that the local fluctuations in energy level dis- 
tribution, as shown in Fig.l for the particular case of SU{1\1) spin chain, get cancelled 
very rapidly whenever we take the cumulative sum of such distribution. Furthermore, 
for sufficiently large values of N, continuous part of the level density distribution in 
the spectrum of supersymmetric HS spin chain satisfies the Gaussian approximation at 
the same high level of accuracy as in the pure bosonic case. It may be noted that, the 
level density of embedded Gaussian orthogonal ensemble (GOE) also follows Gaussian 
distribution at the limit N ^ oo, provided the number of one-particle states tends 
to infinity faster than N [21]. However, in our case of SU{m\ri) HS spin chain, the 
number of one-particle states (i.e., m + n ) remains fixed for all values of N. 

Next we want to study the NNS distribution in the spectrum of SU{m\ri) HS spin 
chain. To eliminate the effect of level density variation in the calculation of NNS 
distribution for the full energy range, it is necessary to apply an unfolding mapping 
to the 'raw' spectrum [22]. This unfolding mapping may be defined by using the 
continuous part of the cumulative level density distribution. We have already seen 
that, for the case of SU{m\n) HS spin chain, the continuous part of cumulative level 
density is given by G{E) (5.2) with a high degree of approximation. So we transform 
each energy i — I,-- - ,Z, into an unfolded energy = G{Ei). The function 
p{u) is defined as the density of the normalized spacings Ui — {^i+i — ^i)/^, where 
A = (^i — — 1) is the mean spacing of the unfolded energy. To get rid of local 
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fluctuations occurring in p{u) , again we study the cumulative NNS distribution given 
by P{u) — p{x)dx, instead of p{u). In this context it may be noted that, NNS 
distributions corresponding to the cases of classical GOE as well as embedded GOE 
obey the Wigner's law [23] : 

p{u) — — wexp(— 7rM^/4) . 

On the other hand, from the conjecture of Berry and Tabor one may expect that 
the NNS distribution for an integrable model will obey Poisson's law given by p{u) = 
exp{—u) [24]. However, it has been found that the NNS distribution for SU{m) bosonic 
HS spin chain does not follow either Wigner's law or Poisson's law within a wide range 
of [9]. Instead, the cumulative NNS distribution for this bosonic spin chain can be 
well approximated by a function like 

='t;"[l -7(1 --y)^], (5.12) 

where v — u/umax with Umax being the largest normalized spacing, a and {3 are two free 
parameters taking values within the range < a, < 1, and the value of 7 is fixed by 
requiring that the average normalized spacing be equal to 1. In our study we also find 
that, NNS distribution in the spectrum of SU {m\n) HS spin chain does not follow either 
Wigner's law or Poisson's law within a wide range of A^. In particular, it is observed 
that the slope of cumulative NNS distribution diverges for both -u — and u u^ax, 
which can not be explained from Wigner's or Poisson's distribution. Furthermore, we 
find that the cumulative NNS distribution for SU{m\n) HS spin chain can be fitted 
well by P{u) in cqn.(5.12) within a range of N. For example, in the particular case 
of SU{1\1) spin chain, it is checked that P{u) agrees well with P{u) within the range 

< 20. In Fig.4, we plot such P{u) and P{u) for = 17 lattice sites {umax = 2.626 
in this case) and found a good agreement with MSE = 0.0234 when the values of free 
parameters are taken as a = 0.39 and j3 = 0.29. However, it is possible that the 
appearance of such non-Poissonian NNS distribution in the spectrum of SU{m\n) HS 
spin chain is an artifact of finite-size effect, which requires further investigation. 

Finally we want to make a comment about the behaviour of parameters fi in 
eqn.(5.8) and a in eqn.(5.11) under the exchange of bosonic and fermionic spin de- 
grees of freedom. Since s s and t ^ —t under this exchange, we find that a remains 
invariant and n changes to p, given hy p, = ^^^^ ~ A* • It is interesting to note that 
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this relation between /i and /i can also be obtained by applying eqn.(4.12): 




N{N^ - 1) 
6 



-II. 



This agreement clearly gives a support to our conjecture (4.11). By using this conjec- 
ture we have found in Sec. 4 that, the spectrum of SU{n\m) spin chain can be obtained 
from that of SU{m.\n) spin chain through an inversion and overall shift of all energy 
levels. Since none of these operations change the standard deviation of level density 
distribution, a should take the same value for SU{'m\n) and SU{n\'m) HS spin chain. 
Hence, the observation that a in eqn.(5.11) remains invariant under the exchange of 
bosonic and fermionic spin degrees of freedom, is also consistent with our conjecture 



Here we derive an exact expression for the partition function of 5'C/(m|n) supersym- 
metric HS spin chain by using the freezing trick and also study some properties of 
the related spectrum. For applying the freezing trick, we consider a spin CS model 
containing an anyon like representation of the permutation algebra as spin dependent 
interaction. We find out the complete spectrum of such spin CS model including the 
degeneracy factors of all energy eigenvalues. At the strong coupling limit, this spin 
CS model reduces to the sum of spinless CS model with only dynamical degrees of 
freedom and SU{m\n) supersymmetric HS spin chain. Consequently, by factoring out 
the contribution due to dynamical degrees of freedom from partition function of this 
spin CS model, we obtain the partition function of SU{m\n) supersymmetric HS spin 
chain. By using this partition function, we study the motif representation for SU{'m\n) 
HS spin chain and find that, due to the lifting of a selection rule, some additional en- 
ergy levels appear in the spectrum in comparison with the case of SU (m) bosonic spin 
chain. 

By using Mathematica we observe that, the partition function of SU{m\n) HS 
model satisfies the duality relation (4.11) for many values of m, n and A^. This obser- 
vation motivates us to conjecture that this duality relation, involving the interchange 
of bosonic and fermionic spin degrees of freedom, is valid for all possible values of m, 
n and N . It would be interesting if this duality relation can be proved analytically 



(4.11). 



6 Conclusion 
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by using the motif representations and skew- Young diagrammes associated with the 
Y{gl[m\n)) quantum group. Furthermore, it is known that, the partition functions of 
SU{m) and SU{m\n) Polychronakos spin chains are intimately connected with Rogers- 
Szego (RS) polynomial [8,15], which appears in the theory of partitions [25]. Since, HS 
spin chain share the same quantum group symmetry with Polychronakos spin chain, it 
might be promising to investigate mathematical structures connected with the parti- 
tion functions of SU{m) as well as SU{m\n) HS spin chain and explore whether some 
new RS type polynomials can be generated in this way. 

By using the partition function of SU{m\n) HS spin chain, we study global prop- 
erties of its spectrum like level density distribution and NNS distribution. It is found 
that, similar to the case of SU{m) bosonic HS spin chain, continuous part of the en- 
ergy level density satisfies the Gaussian distribution with a high degree of accuracy 
for sufficiently large values of N. We also derive exact expressions for the mean value 
and the standard deviation which characterize such Gaussian distribution. It would be 
interesting to provide an explanation for this behaviour of energy level density distri- 
bution in the framework of random matrix theory and explore whether the underlying 
quantum group symmetry of HS spin chain plays some role in this matter. 
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Appendix A. Evaluation of trace formuleis 



Here we shall derive the trace formulas (5.4a,b,c,d) by assuming that |q;i . . . a^) (with 
aj e [1, . . . , s]) are orthonormal set of vectors. Since the trace of identity operator 
is given by the dimension of the Hilbert space, eqn.(5.4a) is really a trivial relation. 
Using eqn.(2.6), it can be shown that 

{ai . . .ai . . .aj . . .aN\r'>j '\ai . . . ai . . . aj . . . aN) = (-1)^ 'Oa.aj, 

where e{ai) = (1) when ctj is a bosonic (fermionic) spin. With the help of above 
equation, we derive the trace relation (5.4b) as 



tr 



(m|n) 



(ai . . . CKj . . . ctj . . . ajvl-P/j"'**^ |q;i . . . ai . . . aj . . . a^) 

s g 

^(ai,aj) ^ l5Hn)| \ 

= [ai . . .ai . . .aj . . .aNi^ij \ai . . . ai . . . aj . . . aN/ 

ai...ajv=l oti=l 
s 

E{ai ,aj) 
represents summation over all spin components cci, . . . , a^v 

except ai and aj. 

Next, by using eqn.(2.6), it is found that 

{ai...ai...aj...ai... a^l^'lj ' F>i '\ai . . .ai . . .a^ . . .ai . . . aN) = Oa^ajOa.ar 

Applying the above equation, we obtain a trace relation in eqn.(5.4c) as 

tr _p.(™l")_p(™l") 

ij il 
s 

= {ai . . .ai . . .aj . . .ai . . . aN\r^ij -r.) '\ai . . . ai . . . aj . . . ai . . . aN) 



* s 

S . 



Other trace relations in (5.4c) can be proved in a similar way. 
By using eqn.(2.6), it is also found that 

[ai . . . ai . . . aj . . . ak . . . ai . . . aN] ^ij -r^i \0(i . . . ai . . . aj . . . a^ . . . ai . . . aN) 



('_l^e(a<)+eK) A A 
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With the help of this equation, we obtain the trace relation (5.4d) as 

tr l^p^'"!")^^™!") 

s 

[Oil . . . ai . . . aj . . . ■ ■ ■ Oil ■ ■ ■ ^Nl-r^ij -r^ki l*^! ■ ■ ■ . . . aj . . . . . . ai . . . a^) 

* s 
ai...ajv=l Q:i,afe=l 

Appendix B. Evaluation of summation formulas 

Here we briefly describe the way of calculating known summation formulas (5.7a) and 
(5.7b) [9], and subsequently present our derivation for new ones like (5.7c) and (5.7d). 
From the work of Calogero et. al [20], it is known that 

> — ir-^ = , B-1 

and 

V 1 _ (iv^-i)(iv^ + ii) 

^sin^(f) 45 ■ ^ ^ 

Using the translational invariance on a circular lattice and summation relation (B-1), 
one can obtain eqn.(5.7a) as 



N-l 



1 ^ 1 _N 1 _ N{N^ - 1) 



Similarly, by using (B-2), one obtains eqn.(5.7b) as 



1 ^ 1 _ ^ 1 _ A^(A^^-l)(Ar2 + ii) 

'"2|-sin^(e.-e,) "^^^h?M~ 90 ■ 

For the purpose of calculating R2 in eqn.(5.5c), we note that Rq in eqn.(5.5a) can 
be expressed as 



1 

k^i sm\Ck-Ci) ' ^ sin^(Ci - Cr) ' ^ sin^(^j - ^^) sm^{Ci-Q' 



Substituting the value of Rq given in (5.7a) to the above relation and also using (B-1), 
we find that 



L 



(Ar2 - l)(7V-4) ^ 1 



sm\^k-^i) 6 ' sin2(ei-0 

{k,lj^i,j) 
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By substituting this expression to R2 in eqn.(5.5c), and subsequently using eqns.(5.7a) 
as well as (5.7b), we derive the value of R2 given in eqn.(5.7c) as 



R2 



V- 

^-^ CI 



• sin (6 - 



(iV2 - l)(iV-4) 
6 



+ 



1 



sin i^i-Q 



N{N^-lf{N-4:) NiN"^ -1){N^ + 11) 



36 



+ 



90 



Finally, by substituting the values of Ro (5.7a), Ri (5.7b) and R2 (5.7c) to the relation 
(5.6), we easily obtain the value of Rs appearing in eqn.(5.7d). 
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Figure 1: Energy levels and degeneracies T)^^^^\Ei) of the SU{1\1) HS spin chain 
for N = 15. 
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Figure 2: Continuous curve represents the error function G(E) and crosses represent 
the cumulative distribution function F(E) (at its discontinuity points) for 5^7(111) spin 
chain with = 15. 
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Figure 3: Left continuous curve represents G(E) for S'?7(l|2) spin chain with N — 15 
and the corresponding F(E) is given by crosses. Middle continuous curve represents 
G(E) for 5'6'^(2|1) spin chain with = 15 and the corresponding F(E) is given by dots. 
Right continuous curve represents G(E) for SU{3) spin chain with N — 15 and the 
corresponding F(E) is given by triangles. 
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Figure 4: Dotted curve represents the cumulative NNS distribution P{u) for SU 
spin chain with = 17 and the continuous curve represents the approximate distribu- 
tion function P{u). 
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